typeset using JPSJ.sty <ver.0.8> 



Critical Dimerization Strength of the Quarter-Filled t—J Model 



00 



(N 



I 

o 



> 

m 
in 
o 

00 



o 
o 



X 



Satoshi NiSHiMOTO^ and Yukinori Ohta^'^ 

^ Graduate School of Science and Technology, Chiba University, Inage-ku, Chiba 263-8522 
^Department of Physics, Chiba University, Inage-ku, Chiba 263-8522 

(Received May 1, 1998) 

An exact-diagonalization technique on small clusters is used to study the ground state of the 
dimerized t—J model at quarter filling. The equal-time charge and spin correlations, charge 
and spin gaps, and binding energy are calculated for the two-dimensional lattice with a spatial 
dimerization pattern corresponding to organic conductors k-(BEDT-TTF)2X. We show that 
competition between the effects of two types of dimerization, i.e., dimerization in the hopping 
integral and dimerization in the exchange interaction, leads to the insulator-superconductor 
transition in the ground state of the model. The phase diagram is thereby presented. 
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It has been emphasizedB* that quasi-two-dimensional 
(2D) organic conductors k- (BEDT-TTF) 2X have some 
characteristics similar to cuprate superconductors; apart 
from 2D nature of electron conduction, strong elec- 
tron correlations, and Mottjdnsiilator to superconductor 
transition|^r|ecfint ^'^C-NMRa'ErQ' and specific-heat mea- 
surementsO^ElEP have suggested that the superconduct- 
ing phase is anisotropic with d-wave-like nodes in the 
gap function and the behavior of (TiT)^^ is somewhat 
similar to the spirbjgap phenomena observed in cuprate 
superconductors Also noticeable is the very short co- 
herence length of ^o=25^30AJ13' which is comparable to 
that for cuprate superconductors. 

Motivated by these intriguing experimental data on 
the materials and some theacetical calculations on the 
dimerized Hubbard models JIj^' we here study the ef- 
fects of dimerization on the ground state of the quarter- 
filled t—J model, another representative model for 
strong electron correlations, by focusing on its insulator- 
superconductor transition. We will show that due 
to competition between two types of dimerization the 
ground-state phase of the model is divided into a Mott- 
insulating phase and a singlet-pairing phase; i.e., there is 
a critical dimerization strength. Although applicability 
of the model to the organic compounds is not yet fully 
worked out, we believe that the results shown here will 
help one understand a possible mechanism of supercon- 
ductivity in strongly correlated electron systems such as 
K-(BEDT-TTF)2X. 

The dimerized t—J model we study is defined by the 
Hamiltonian 

H=-t,Y, (cLc.a + H.c.) -t^Y, {i^^^ha + H.c.) 



) + J2^(Sfc•S^- 



r^fer^^ 



)(0.1) 



where cl^—cl^i^^n-i^^) is the constrained electron- 
creation operator at site i and spin a {—"l, j), Si is the 
spin-i operator, and Ui is the electron- number opera- 



tor; we refer to the fermionic particle as 'electron', which 
corresponds to the hole in the real organic compounds. 
(ij) stands for nearest-neighbor bonds with parameters 
ti and Ji and (kl) for those with parameters t2 (>ii) 
and J2 {>Ji)- We assume the lattice of a dimerization 
pattern simulating the 2D conducting plane of k-(BEDT- 
TTF)2X compounds (see Fig. 1(a)). Note that the model 
tends to the usual square-lattice t—J model when there 
is no dimerization (ij=ia and Ji=J2), of which much 
study has been made,Ej'lHl' whereas in the limit of strong 
dimerization, the model represents an assembly of iso- 
lated dimers. We retain the relations between parame- 
ters t and J obtained from perturbation, i.e., Ji=4i^/J7 
and J2=4t|/t/, in order to reduce the number of pa- 
rameters, where U is the corresponding on-site Hubbard 
interaction. We thereby keep a relation Ji/ J2={ti/t2)^ ■ 
We thus have three independent parameters, and if we 
take ti as a unit of energy, then we have two, for which 
we will take parameters representing t-dimerization and 
J-dimerization (of which a specific definition is given be- 
low). Here we restrict ourselves to the case of quarter fill- 
ing to simulate the situation of the organic compounds. 

We introduce two types of dimerization: one is 
the dimerization of hopping integral, which we call t- 
dimerization, and the other is the dimerization of ex- 
change interaction, which we call J-dimerization. We 
define a parameter 

t2~ti 



td 



(0.2) 



for the strength of i-dimerization and a parameter 

J2-J1 



Jd 



(0.3) 



for the strength of J-dimerization. We also take J2/t2 
as a measure of the strength of J-dimerization because 
if we keep id constant (> 0) then J2/i2 esc Jd- We note 
that the t-dimerization has the effect leading to a repul- 
sive interaction among electrons that acts when differ- 
ent spins (1" and [) come in a single dimer and that the 
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J-dimcrization has the effect promoting the spin-singiet 
formation between spins coming in a single dimer. Thus, 
the competition between the effects of these two dimer- 
ization may result in the following situation as illustrated 
in Fig. 1(b); when i-dimerization is dominant, the system 
can be a Mott insulator (because we have an effective 
half- filled band with the 'on-site' repulsion UdimcJ^), 
when J-dimerization is dominant, the system can be a 
singlet-pairing state, and when J2/t2 is very large, the 
system will be phase separated. 



Fig. 1. (a) A 2D lattice structure of k-(BEDT-TTF)2X. Tiie 
nearest-neighbor bonds have either ti and Ji (thin solid hne) 
or t2 and J2 (bold line). The unit cells, each of which contains 
four sites, are indicated by dashed lines, (b) Schematic ID rep- 
resentation of the electronic states: from top to bottom, Mott 
insulator, singlet-pairing state, and phase separation are illus- 
trated. 



In the following, we will present numerical evidences 
that the competition indeed leads to the insulator- 
superconductor transition in the ground-state phase dia- 
gram of the dimerized t—J model. We employ a numeri- 
cal exact-diagonalization technique on small clusters. We 
use clusters of the size 1x4 and 2x2 unit cells (16 sites) 
with periodic boundary conditions (see Fig. 1 (a)). We 
find that nearly the same results are obtained for these 
two clusters, so that hereafter we will show the results 
only for the cluster of 2x2 unit cells. 

Let us first discuss the equal-time charge and spin 
correlations. In Fig. 2, we show the calculated results 
for the charge correlation (riinj) and spin correlation 
(SfSj), where (• ■ •) denotes the ground-state expecta- 
tion value. For comparison, we also show the results 
for the t—J cluster without dimerization {td=Jd=0) in 
the topmost panels of Fig. 2. It is first of all evident in 
the figures that the effect of dimerizations indeed bring 
strong impact on the intra-dimer charge and spin corre- 
lations. We find that, in the small J2/i2 region where J- 
dimerization is weak, the intra-dimer charge correlation 
becomes closer to zero, indicating a tendency that only 
one electron presents in a dimer. This suggests that the 
effective repulsive interaction presents in a dimer. With 
increasing i-dimerization (td) this tendency appears to 
be much enhanced as clearly seen in Fig. 2. On the other 
hand, in the large J2/t2 region where the J-dimerization 
is large, the intra-dimer charge correlation is approach- 



Fig. 2. Charge correlation (riinj) (left panel) and spin correlation 
{S^Sj) (right panel) for the dimerized t—J model. Intra-dimer 
correlations are indicated by full squares. Topmost two panels 
are for the cases without dimerization. 



ing the value 1 /4 and the intra-dimer spin correlation is 
decreasing to —1/16 as clearly seen in Fig. 2. This result 
indicates that two electrons with opposite spin come in a 
single dimer when J-dimerization is large. Then, in the 
region of intermediate strength of J2/t2, there appears a 
state where the charge correlations among different sites 
have nearly the same value ~0.5/4 and thus there are 
rather small spatial correlations. We note however that 
even in this region the nearest-neighbor spin correlations 
are still significantly antiferromagnetic. 

Next let us examine the single-particle gap and charge 
gap. For the 16-site cluster at quarter filling the single- 
particle gap may be defined by 

Aei = i[[^Gs(5,4)-i;Gs(4,4)] 

- [^Gs(4,4)- J;gs(3,4)]] (0.4) 

where i?GS (-^T i ) i^ '^^^ ground-state energy for the 
cluster of N-^ up-spin and iVj^ down-spin electrons. The 
calculated results are shown in Fig. 3 (a). When there 
is no dimerization, one might expect that the system is 
metallic, i.e., Aci=0 in the entire region of J2/i2 (except 
in the region of phase separation) . We see however that 
the obtained value for Ad is finite and it increases with 
increasing of J2/i2- The reason is, apart from an obvi- 
ous finite-size effect of small clusters, that Ad reflects 
the effect of electron pairing, even in the charge-gapless 
region (as discussed below) . When there is dimerization, 
Ad is finite even at J2 7^2=0 and rapidly increases with 
increasing the strength of t-dimerization. We see that 
at constant td and with increasing J-dimerization J2 /t2 , 
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Fig. 3. (a) Single-particle gap Aci 7*2, (b) charge gap Ac2 7*2, (c) 
spin gap As 7*2) and (d) binding energy ^.-gji^ as a function of 



the gap first decreases, reaches a minimum, and then in- 
creases, i.e., there appears a minimum value in Ad. The 
value of Jilt2 at which this minimum occurs is defined 
as {J2lt2)c- The minimum position (J2/t2)c increases 
with increasing t-dimcrization. 
The charge gap defined by 

n[^Gs(5,5)-£;Gs(4,4)] 



Ac2 = 



[i?Gs(4,4)-i?Gs(3,3)] (0.5) 



may also be introduced. The calculated results are shown 
in Fig. 3 (b). When there is no dimerization, we find the 
result Ac2<0, contrary to a simple expectation. This 
may however be understood if we assume that the attrac- 
tive interaction between pairs (or some tendency toward 
phase separation), of which effect Ac2 reflects, exists in 
the singlet-pairing region. When there is dimerization, 



the charge gap Ac2 opens for td>Q and increases rapidly 
with increasing td- With increasing J2/i2, the gap de- 
creases first and at some J2/i2 value the rate of the de- 
crease changes, and at the same time Ac2 changes sign 
where the gap closes. We may write this J2/<2 value 
as (J2 7^2)0 because the value agrees with the value of 
{J2li'i)c defined by using A^i. Above (^2/^2)^ the gap 
decreases as 321^2 increases, which may again be due to 
the effect of attractive interaction between pairs. 

Let us now consider the spin gap, which is defined by 



A, = Sgs(5,3)-£;gs(4,4). 



(0.6) 



The calculated results are shown in Fig. 3 (c). We find 
that the spin gap opens at J2/t2'^( ^2/^2)0 and increases 
rapidly with increasing J2/t2 (if we assume that a small 
As>0 at J2/t2<{J2/t2)c is a finite-size effect). We also 
find that with increasing td the size of the gap decreases 
and the value {J2/t2)c increases. Here we should em- 
phasize that the (^2/^2)0 values estimated from the be- 
haviors of Aci, Ac2, and A<, are all consistent and have 
the same td dependence. Thus it seems quite reasonable 
to assume that the spin gap Aj is finite in the entire 
charge-gapless region, where the electrons form a bound 
state as shown below. 

The binding energy may be estimated by 

Ab = [£;gs(5,5)-Sgs(4,4)] 

-2[£;gs(5,4)-£;gs(4,4)] (0.7) 

where As is negative if the state of two electrons mini- 
mize its energy by forming a bound state. A possibility 
of superconductivity may then be indicated. The calcu- 
lated result for Ab is shown in Fig. 3 (d). We find that 
|Ab| is smaU for J2lt2<{J2/t2)c but starts to increase at 
around (^2/^2)0, which is consistent with the expectation 
that the bound state of electrons is formed in the region 
of charge-gapless and spin-gapful phase. We also note 
that, with increasing i-dimerization, the binding energy 
is suppressed and the value of {J2/t2)c increases accord- 
ingly. We should however be careful because when the 
value of ^2/^2 is too large, the system will be phase sep- 
arated. We estimate the parameter region of phase sep- 
aration by examiiiiag if the value of the compressibility 
becomes negativcEJ* It is known that the phase separa- 
tion actually occurs-at around 72/^2—2.5—3 when there 
is no dimerization.119 We thus examine its dependence 
on the t- and J-dimerizations; the obtained results are 
summarized in the phase diagram shown below. 

Based on the calculated results given above, we ob- 
tain the phase diagram of the model on the parameter 
space of the two dimerizations. The result is shown in 
Fig. 4 where the horizontal axis is the t-dimerization td 
and vertical axis is J2 /t2 , the measure of J-dimerization. 
The value of Jd is also plotted. We see that when the 
i-dimerization is predominant over the J-dimerization, 
the system is a Mott insulator (the effective half-filled 
band with iJdimor), whereas when the J-dimerization 
is predominant over the t-dimerization, the syjStem is 
in a singlet-pairing (or superconducting) statctiP The 
boundary between these two phases approaches J2/i2=2 
in the limit of strong t-dimerization. This phase diagram 
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Fig. 4. Ground-state phase diagram of the 2D dimerized t—J 
model at quarter filling obtained for the 2x2 unit-cell cluster. 
Values of are indicated as a contour (dashed line). We put 
some values of the spin gap As/t2 in the singlet-pairing region. 



is thus quite consistent with the expectation discussed in 
the beginning of the paper and illustrated in Fig. 1 (b). 

It should be noted here that the phase diagram is ob- 
tained without finite-size scaling analysis of small-cluster 
data, as it is not feasible in any 2D-model calculations. 
We point out however that almost the same results are 
obtained for the 2x2 and 1x4 unit-cell clusters, that 
indications for the existence of the critical dimerization 
strength are robust, and that the proposed physical pic- 
ture is quite obvious. We also note that, in the ID model 
where a finite-size scaling can be done, we have confirmed 
the same mechaiiism for the insulator-metal transition 
actually worksO* We therefore believe that the numer- 
ical results presented in this paper reflects the reality of 
the model at least qualitatively. 

Finally, let us examine the recent experimental 
data and some other calculations by taking k-(BEDT- 
TTF)2X as an example and see if we can find 
any relevance to our model system. In k-(BEDT- 
TTF)2Cu(NCS)2, the hopping integrals ti and t2 for 
holes are reported to be ti~Q.l eV an(i-£2— 0.25 eV from 
the extended Hiickel band calculation,ll3) which leads to 
the strength of <-dimerization td~1.5. The Coulomb in- 
teraction U on a BEDT-TTF molecule is reported to 
be ~1.0 eV, which is deduced from the analysis of the 
pressure dependence of the Kaight shift based on the 
random-phase approximation.lij' Thus the Coulomb in- 
teraction may be 4—10 times larger than the hopping 
integrals. Then if we assume that the exchange inter- 
actions can be estimated by the perturbation expression 
J=U'^/U we have Ji/ti~0.4 and J2/t2^1.0. It follows 
then that the obtained point {id, J2/t2) falls on the Mott- 
insulating region of the phase diagram of Fig. 4. Now it 
is interesting to note that the effect of pressure is to 



decrease ^2/^1 (and simultaneously increase J2/J1); an 
optimistic expectation would thus be that the system 
which is located in the Mott-insulating phase at ambi- 
ent pressure shifts to the upper-right direction in Fig. 4 
and can be in the singlet-pairing (or superconducting) 
phase with increasing pressure, as it actually is in ex- 
periment. It is also interesting to note that the critical 
dimerization strength (J2/i2)c— 1-7 obtained at td—l.^ 
(see Fig.|_4) is roughly consistent with the Hartree-Fock 
estimatdiiP of the critical repulsion (?7dimcr/^2)c— 1-0 in 
the dimerized Hubbard model because, if we use the ex- 
pression L'dimci~2i— J of an isolated t—J dimer, we have 
( J2/i2)c— 1-0. Further studies will however be required 
to clarify whether there is any continuity between the 
present dimerized t—J model and more realistic dimer- 
ized Hubbard models, and also whether such large J/t 
values are realized in the organic compounds k-(BEDT- 
TTF)2X. 

In summary, we have examined the 2D dimerized t—J 
model at quarter filling by using an exact-diagonalization 
technique on small clusters and have shown that there is 
a critical dimerization strength which divides the ground- 
state phase of the model into the Mott-insulating phase 
and singlet-pairing phase; either of the two phases is real- 
ized depending on the strength of t- and J-dimerizations. 
We have thereby obtained the ground-state phase dia- 
gram of the model in the parameter space of the two 
types of dimerization. 
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